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We argue that a locally time-reversal symmetry breaking 
state can occur at Josephson junctions between unconven- 
tional superconductors. Order parameters induced by the 
proximity effect can combine with the bulk order parame- 
ter to form such a state. This property is specifically due to 
the intrinsic phase structure of the pairing wave function in 
unconventional superconductors. Experimental consequences 
of this effect in high-temperature superconductors are exam- 
ined. 



A growing number of experiments demonstrates 
that the superconducting state in some of the high- 
temperature superconductors (HTSC) possesses an un- 
conventional order parameter which changes sign under 
a 90°-rotation in the copper-oxide plane 0] . For a tetrag- 
onal system this state is commonly refered to as a d- 
wave state, where the Cooper pairing occurs in a channel 
belonging to an irreducible representation of the point 
group D^fi, which is different from the trivial representa- 
tions, i.e. "s-wave". Many of these experiments use the 
Josephson effect, which allows a direct probe of the phase 
properties of a superconducting order parameter. While 
the Josephson effect in conventional (s-wave) supercon- 
ductors is described only by the phase difference between 
the order parameters on both sides of an interface, the 
situation is more complicated if the pair wave function 
has an intrinsic phase structure like in the d-wave state. 
The sign change of the d-wave state (generic pair wave 
function: ipO^) — cosfc^. — cos ky) between the x- and y- 
direction corresponds to a phase difference of tt. Due to 
this intrinsic phase structure, the relative orientation of 
the pair wave function of two connected superconductors 
and the geometry of the connecting interface are impor- 
tant to determine the phase relation between the order 
parameters on both sides [§-1). For d-wave supercon- 
ductors two types of Josephson junctions can occur, one 
where the interface energy is minimized by a phase dif- 
ference if = (standard 0-junction) and the other with 
(p = TT (TT-j unction) In multiply connected supercon- 
ducting systems the latter type can lead to frustration 
effects and twists of the order parameter phase . Re- 
sulting phenomena like the occurrence of spontaneous su- 
percurrents and the modification of standard interference 
patterns (SQUID) have been exploited in experiments to 
determine the symmetry of the order parameter Q . 



In this letter we add a new feature to the properties of 
interfaces between unconventional superconductors. We 
examine the effect of the reflection of Cooper pairs at 
the interface and the mutual proximity of the two su- 
perconductors, both of which can introduce Cooper pair 
amplitudes in channels of symmetries different from that 
of the bulk pairing state. These additional pairing com- 
ponents can exist only in the immediate vicinity of the in- 
terface and decay exponentially towards the bulk. How- 
ever, their presence has important consequences for the 
interface properties because they modify the Josephson 
effect via their influence on the current-phase relation. 
The admixed order parameters originating from reflec- 
tion and proximity usually prefer a relative phase of 
or TT with the bulk order parameters on both sides of 
the interface, resulting in states conserving time-reversal 
symmetry. We will show here that under certain con- 
ditions the relative phase can change leading to a state 
which breaks time-reversal symmetry at the interface. It 
was shown by Sigrist, Bailey and Laughlin that this 
kind of interface state may explain a recent experiment 
done by Kirtley et al. where fractional vortices have been 
observed at boundaries between differently oriented films 
of the HTSC YBazCuaOT-i (YBCO) §. 

We discuss this problem using the Ginzburg-Landau 
(GL) theory of a d-wave order parameter in a tetrago- 
nal system {D^^h) and consider only one admixed order 
parameter (s-wave) which does not occur in the bulk. 
This choice is motivated by the experimental results in- 
troduced above |0] and the recent finding that the critical 
behavior at the onset of superconductivity suggests a sin- 
gle component bulk order parameter in HTSC ||]. The 
general GL free energy expansion for these two compo- 
nents, "qd (d-wave) and r^s (s-wave), in two dimensions is 
given by 
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with aj{T) — aj{T — Tcj) and ctj, Kj, 71, 72 and K 
are real coefficients (j = d, s). The gauge invariant gradi- 
ent is given by D = V — i(27r/$o)A with A as the vector 
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potential and $o = /ic/2e is the standard flux quantum. 
We consider the properties of an infinite planar interface 
with the geometry shown in Fig.l where the crystal main 
(x-) axis of the side A is pointing perpendicular to the 
interface, while the orientation of the side B is described 
by the angle 9 between the interface normal vector and 
the crystal main (x-) axis (— 7r/2 < < tt/2). In this 
arrangement the effects of interest occur only on side B. 
Thus, we will neglect the detailed description of side A 
and represent it simply by a complex (bulk) order pa- 
rameter rjo- It can have d-wave symmetry, but s-wave 
symmetry would not change any of our conclusions qual- 
itatively. The properties of the interface are described by 
additional interface terms in the GL free energy 



FjF = J,^ dS[gd{e)\Vd\' + m{v*dVs + VdVl) 



(2) 



The first two terms describe the reflection properties of 
the interface (IF) for the side B, i.e.: r]d — > rjd in the first 
and rjd ^ rjs in the second term (we neglect the rjs rjs 
contribution). The latter two terms represent the lowest 
order coupling between the two sides, A and B. By sym- 
metry argument we find that g{9 ± 7r/2) = —9{9) and 
ti{0±Tr/2) = -ti{e). Thus, we choose g = go cos(26') 
and ti = tocos{29) with 50,^0 < as a generic angle 
dependence. We neglect any ^-dependence in t2 because 
none is required by symmetry (^2 < 0). The signs of go, 
to and t2 are chosen by convention defining the specific 
gauge of the order parameter phases. The coefficient gd is 
in general positive and describes the suppression of the 
d-wave order parameter at the interface (see The 
variational minimization scheme for the GL free energy 
includes these interface terms into the boundary condi- 
tions. 

This GL theory contains two competing tendencies for 
the relative phase, cps — (fid, of the order parameters 
r7j(= \r]j\exp{iipj)), {j = s,d). The last three terms in 
Fip (Eq.(2)) are minimized by (/s^ — "ySd = for 16*1 < 7r/4 
and (fs — (fd = for 7r/4 < j^j < tt/2. In either case, the 
state at the interface is invariant under time-reversal op- 
eration. On the other hand, if we assume that 72 > 0, the 
coupling term ?7i^f72 +^71^2^ Eq.(l) favors energetically 
'■Ps~^d = i7''/2, a state which breaks time- reversal sym- 
metry (s-|-id-wave state). (The assumption of positive 
72 is natural, since it would lead to a fully gapped state 
which is energetically more favorable as weak-coupling 
studies of various microscopic models show |^.) In the 
following we study the competition between the tendency 
to keep or to break time-reversal symmetry near the in- 
terface and some physical consequences. 

Although our choice of geometry allows to reduce the 
problem to one spatial dimension {x = x cos 9 + y sin 9 
on side B), these GL equations can only be solved nu- 



merically in general. Fortunately, it is possible to under- 
stand the essential properties of the competiting effects at 
the interface qualitatively by using a simple variational 
ansatz for the order parameter. We assume that as in 
HTSC the transition temperature Ted of the d-wave order 
parameter is finite, while that of the s-wave is very small 
or zero, Tcs <C Ted such that in the bulk only rjd is finite 
for < T < Ted (and 770 on the side A) . In our varia- 
tional treatment r]d shall be real and independent of posi- 
tion on the side B, and 770 of the side A shall have a fixed 
modulus fj and a phase x, rjo = fjexp{ix) {rid{T),fi{T) 
are identical to their bulk values). Hence, this ansatz ne- 
glects the direct effect of the interface on the modulus of 
the d-wave order parameter. The s-wave component 77^ 
is complex and decays exponentially towards the bulk on 
the side B, ris{x) = f]e~^/^ and rj — fj' + if]" ■ The part 
of the free energy (per unit interface area) depending on 
the variational degrees of freedom, ^, 77 and x, is obtained 
straightforwardly 
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'2{g{9)r]dfl' +ti{9)firid cos X (3) 
'hvifj' cosx + fj" sinx)} 



with K± = Ks + a±^^ and a± = Us + (71 ± 72)'?^- To be 
consistent with our assumption that 77s = in the bulk 
we require a± > 0. Further, we choose a_ < a+ (72 > 0). 
By minimizing Fyar we find that ^ is always finite with 
Ks/a+ < < Ks/a-. Some simple algebra leads to the 
following equation for x, 



sinx [cosx + 



{hK+ - 2Cgt2)K^ 
4iiC^72?7% 



(4) 



One set of solutions of this equation is x = and tt 
(sin X = 0) which leads to 77' 7^ and 77" — such that the 
relative phase between 77^ and 77s is or tt. Energetically 
X = (77' > 0) is favored for \9\ < tt/A (0-junction) and 
X = TT (77' < 0) for \9\ > 7r/4 (7r-junction). However, also 
the term in [...] can vanish giving an alternative solution 
of Eq.(4) with x different from the above two limiting 
values. Both 77' and ff are finite in this case. This state is 
two-fold degenerate because the application of the time- 
reversal operation (7/ — > 77*) corresponding to x ~X 
and fj —>■ f]* leads to a different state with the same free 
energy. Consequently, this solution violates time-reversal 
symmetry T. It is easy to see that the energy is indeed 
minimized by this state if 72 is positive. From Eq.(4) it 
is obvious that the T- violating state can only occur if the 
modulus of the second term in [...] is smaller than 1. This 
condition can be satisfied by an appropriate choice of the 
angle 9, because this term is proportional to cos(26') via 
ti and g, which can be arbitrarily small for 9 close to 
±7r/4. 
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From the above consideration the condition for the in- 
stabiUty of the time-reversal invariant state is obtained 
as 



(^0°' and (^q"'' , respectively) , we find a continuous change 

of i-p from (^q"' to ((Sq^"* on a length scale Aj. On dis- 
tances much larger than A,/ from the border of the seg- 
ments <~p is (^q"'"' in segment (j). The flux associated with 
the variation of ip is located in the vicinity of the bor- 
der and is $ = (ip^^^ — (p'^^)^o/2Tr. Obviously $ is not 
an integer (half- integer) multiple of $0i but can have an 
arbitrary value. It cannot be determined by simple topo- 
logical arguments, but depends on the specific properties 
of the interface. Because the phase twist at the border 
between segment (a) and (b) corresponds to a winding of 
the phase ipd by a fractional multiple of 27r, this flux line 
may be considered as a fractional vortex [ p^ . Similarly 
fractional vortices can occur on a homogeneous interface 
due to the twofold degeneracy of the T-violating interface 
state. Two types of domains with ipQ — +x and — x, re- 
spectively, are possible. If present on the same interface, 
they are separated by domain "walls" . There we find a 
kink of the width Aj in ip connecting the two values of 
ipo and yielding a flux $ — ±x^o/'^ (see also flo|]). 

We now propose an experiment to test our picture. Let 
us consider a thin film consisting of two parts with orien- 
tations of their crystalline axes different by 45°. The 
boundary shall be circularly curved (see Fig. 3). The 
point P would correspond to an interface with angle 
9 = tt/A and changes continuously as we turn away from 
P. Hence, along the interface we scan through a certain 
range of angles 9. If time-reversal symmetry were con- 
served everywhere, we would see a vortex with $ = $o/2 
at P, because P separates two interface segments where 
one has (po = and the other = tt |2[[Ti||. The field 
distribution of the vortex becomes more localized if we 
cool the system because the penetration depth Aj along 
the interface shrinks with lowering temperature. On the 
other hand, if our scenario is correct we expect that the 
phase ipo varies continuously from to tt near P within 
a certain range of 0. This range extends with lowering 
temperature (see Fig. 2), such that the resulting field dis- 
tribution becomes more extended further as the system 
is cooled, while the total flux always equals to $o/2. The 
observation of the field distribution by a (magnetic) mi- 
croscope should allow to prove or disprove our scenario 
based on these qualitative properties. 

In summary, we pointed out that important properties 
of an interface between d-wave or other unconventional 
superconductors result from the intrinsic phase structure 
of the pairing wave function V'(k) which cannot be found 
in conventional (s-wave) superconductors. While some 
interface coupling effects favor a T-invariant state, in the 
bulk rather a T-violating state is prefered. In the compe- 
tition between these two tendencies it is important that 
the interface coupling effects can be arbitrarily small by 
choosing the interface geometry appropriately, i.e., the 
crystalline orientation of the two connected superconduc- 
tors. Thus, proximity-induced order parameter compo- 



cos(26') = ± 



th2fliT*)vdiT* 



(5) 



where T* (< Td) denotes the critical temperature of 
the second order transition to the T-violating state for a 
given angle 9. The denominator on the right hand side 
of this equation is only weakly dependent on T* close to 
Ted so that T* is large for near ±7r/4, and is Ted — 
±7r/4. We show a schematic phase diagram, T versus 9, 
corresponding to Eq.(5) (Fig. 2). The T-violating phase 
appears in a certain range of around 7r/4 and the region 
of this phase becomes wider with decreasing temperature. 
The specific shape of the phase boundary line depends on 
properties of the system, in particular, of the interface. 

Let us now investigate important consequences of the 
existence of a T-violating interface state, which result 
from the properties of its Josephson current-phase rela- 
tion. 



(6) 



where ip is the phase difference between the order param- 
eters rjQ and r]d and .g((p) is a 27r-periodic function of p 
{\9{^)\ ^ 1) The specific form of g is not important 
here and will be discussed elsewhere. The interface state 
with minimal energy corresponds to / = and its twofold 
degeneracy implies that two functions, and (?_, ex- 
ist with g+{x) — and 5-(— x) — 0, each belonging to 
one of the two states. There are two possible situations: 
I) g+{(p) 7^ g-if) describing both different (meta)stable 
states as a function of and II) g+{p) — g-{^p) where 
the two T-violating states can be adiabatically connected 
with each other by changing Lp from -1-^ to — The spa- 
tial variation of (p along the interface (with coordinate 
denoted as x') is described by a generalized Ferrel-Prange 
equation which is identical to a Sine-Gordon equation for 
standard Josephson junctions {g(pi) — sintp), 



A-2 



(7) 



with A J = (<I>oc/87r^(i/c)^/^ as the characteristic length 
scale (d: the effective magnetic width of the interface 
[^). As in conventional interfaces the spatial varia- 
tion of Lp yields a finite local magnetic flux (/)(x') = 
{^o/2TT)dip{x') /dx' . The integrated magnetic flux be- 
tween two points, x'2 and x'l, on the interface is given by 
$ = (^(4)-(p(x;))$o/27r.^ 

As a consequence, a spatial variation of the interface 
properties can lead to the occurrence of magnetic field 
distributions. In particular, at the border between two 
different homogeneous interface segments, (a) and (b), 
characterized by their specific values of (po (denoted as 
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ncnts can combine with the bulk order parameter to form 
a T-violating interface state. We emphasize that such 
conditions cannot be satisfied in conventional (s-wave) 
superconductors, in general. Therefore the proximity- 
induced time-reversal breaking interface state is a con- 
sequence of unconventional superconductivity. For the 
sake of simplicity we restricted ourselves to the case of a 
system with tetragonal crystal field symmetry. We note 
here, however, that a slight orthorhombic distortion of 
the lattice (as it is present in real materials) does not 
change our conclusions qualitatively. It leads to some 
minor modifications which are beyond the scope of this 
letter and will be discussed elsewhere in detail. 
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FIG. 2. Schematice phase diagram in terms of temperature 
and crystal orientation of side B. The shaded region denotes 

the T-violating phase separating the 0- and the yr-junction 
phases. Within the shaded region x changes continously. 

FIG. 3. Schematic arrangement of superconducting films 

for a test-experiment. The direction of shading indicates the 
crystalline orientation. The region A may be a d- or an s-wave 
superconductor, while B must have d-wave symmetry. The 
point P corresponds to the angle ^ = 7r/4 and 9 = ■k/A — 9. 
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FIG. 1. Interface between d-wave superconductor A and 
B. The shading on either side indicates the direction of the 
crytalline a;-axis (the z-axis points out of the plane). 
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